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Abstract 

R measurement data taken with the BESII detector at center-of-mass energies 
between 3.7 and 5.0 GeV is fitted to determine resonance parameters (mass, total 
width, electron width) of the high mass charmonium states, ^(3770), 0(4040), 
■0(4160) and -0(4415). Various effects, including the relative phases between the 
resonances, interferences, the energy-dependence of the full widths, and the initial 
state radiative correction, are examined. The results are compared to previous 
studies. 

1 Introduction 

The total cross section for hadron production in e + e~ annihilation is usually parameterized in 
terms of the ratio R, which is defined as R = a(e + e~ — > hadrons)/<r(e + e~ — ► /i + /i _ ), where the 
denominator is the lowest-order QED cross section, a(e + e~ — ► fi + yu~) = crjL = Ana 2 /3s. The 
measured R values are consistent with the three-color quark model predictions. At the open 
flavor thresholds where resonance structures show up, measurements of R value are used to 
determine resonance parameters. For the high mass charmonium resonances, the 0(3770) was 
measured by MARK-I PQ, DELCO [2], MARK-II [3] and BES @][5]; the 0(4040) and 0(4160) 
were measured by DASP (Bj; and the -0(4415) was measured by DASP [6] and MARK-I (7J. 
There were also some other measurements of R values as reported in Refs. [U [91 [10], but 
no attempt was made to determine resonance parameters. The resonance parameters in the 
Particle Data Group (PDG)'s compilation remained unchanged for more than 20 years up 
to the 2004 edition [IT]. The resonance parameters for the three high mass resonances were 
updated by PDG2006 [J2], based on K. Seth's evaluation [TO] using combined BESII [U] and 
Crystal Ball [TO] data. 
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The most recent and precise R measurements between 2-5 GeV were made by BESII |14j . 
Experimentally, R for both the continuum and the wide resonance region is given by 



R, 



exp 



where Nj^ s d is the number of observed hadronic events, Nf, g is the number of the residual 
background events, L is the integrated luminosity, (1 + 5 b s ) is the effective correction factor 
of the initial state radiation (ISR) [10] [15] , e^ad is the detection efficiency for hadronic events 
determined by the Monte Carlo simulation without bremsstrahlung being simulated, and €t rg is 
the trigger efficiency. The determination of R values and resonance parameters are intertwined; 
the factor (1 + 5 &s) in Eq. ([T]) contains contributions from the resonances and depends on the 
resonance parameters. Therefore, the procedure to calculate (1 + 5 b s ) requires a number of 
iterations before stable results can be obtained. 

In this work, we perform a global fit over the entire center-of-mass energy region from 3.7 
to 5.0 GeV covering the four resonances, "0(3770), 0(4040), ^(4160) and 0(4415), and include 
interference effects among the resonances. We also adopt energy-dependent full widths, and 
introduce relative phases between the resonances. Finally, the new R values due to the updated 
resonant parameters are compared to the published ones [14] . 

2 Fitting of the resonant parameters 

The phenomenological models and formulas used in our fitting are briefly described below. 

2.1 Breit-Wigner form 

The relativistic Breit-Wigner amplitude for e + e~ — ► resonance — ► final state / is 



T/(\V) V r l) 



W 2 - AG + iM r T. 



-e 



where W = E cm = y/s is the center-of-mass energy, the index r represents the resonance to 
be studied, M r is the nominal mass, r r is the full width, is the electron width, r£ is the 
hadronic width for the decaying channel /, and 8 r is the phase. 

The high mass charmonia decay into several two-body final states /. According to the 
Eichten model [16] and experimental data [T7], the decay channels (including their conjugate 
states) are: 

0(3770) DD; 

0(4040) => DD,D*D*,DD*,D S D S ; 

0(4140) DD,D*D*,DD*,D S D S ,D S D*; 

0(4415) DD,D*D*,DD*,D S D S ,D S D* S ,D*D*,DD 1 ,DD%. 



3 



The total squared inclusive amplitude of the resonances is the incoherent sum over all 
different decay channels /, 

i^i 2 = EiE r /wi 2 - ( 3 ) 

/ r 

The resonant cross section expressed as the R value is then given by 

p _ °>es _ 127T . 2 f a\ 

r Wes — n — \ 1 res\ ■ V / 

°W S 

2.2 Energy-dependence of the full width 

The full width of a broad resonance depends on the energy. A phenomenological model derived 
from quantum mechanics is used to describe the behavior of Tf,(W), which depends on the 
momentum and the orbital angular momentum L of the decaying final state [18] . 

where, T r is a parameter to be determined by fitting experimental data, Zf = pPf, p is the 
radius of the interaction with the order of a few fermis (the value is insensitive to the physical 
results), Pf is the decay momentum, and 2M r /(M r + W) is a relativistic correction factor. 
The energy-dependent partial wave functions P>l are given in Ref. |18j : 

B = l, B 1 = l + Z 2 , B 2 = 9 + 3Z 2 + Z 4 , 

B 3 = 225 + 45Z 2 + 6Z 4 + Z 6 . (6) 

When the resonance decays to several hadronic channels, the total hadronic width is the 
sum of all its partial widths, 

T h r ad {W) = Y J ?l(W). (7) 

/ 

The total width of the resonance r is expressed as 

T r (W) = 3T e r e + T^ ad (W), (8) 

where the universality of the leptons is used, i.e. T^ e = rjf = T^ T , and T^ e is the experimental 
electron width, which includes the contribution from the vacuum polarization effect. 



2.3 Continuum background 

The contribution of continuum background originating from initial light quark pairs uu, dd 
and ss, which is written as Rqcd> ma y ^ e predicted by pQCD above 2 GeV|12|. Being close 
to the production threshold, the continuum open-cc background can only be described by 
phenomenological models or empirical expressions. Since there are many possible channels 
above the open-charm threshold, and their cross sections are expected to vary smoothly, we 
parametrize, for simplicity, the continuum charm background with a second order polynomial, 

flW = C + Ci ( W - 2M D± ) + C 2 ( W - 2M D± ) 2 , (9) 

where Co, C\ and C2 are free parameters, and M D ± is the mass of the lightest meson D . 
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2.4 Fitting scheme 

We fit the experimental data with MINUIT [19] using a least squares method, with x 2 defined 



as 

x2 = ^ [fcRexp(W i )-R the (W i )] 2 | (/c-1) 2 ^ (1Q) 
i [fc^Rexp] 2 °" c 

where Wi stands for the energy of the measured point. The experimental quantity 

R exp = N P ~ Nb9 , fll) 



and the corresponding theoretical quantity 

Rthe = (1 + 5 bs)Rthe, (12) 

where, 

Rthe — Rcon ~i~ Rresi (1^) 

and 

Rcon = R(QCD) Rcon- (1^) 

ARexp in Eq. ()10p is the combined statistical and uncommon systematic errors of R e xp(Wi)\ 
the error common to all the points a c (~ 3.3%) is not included. In each iteration, the resonant 
parameters in the calculation of (1 + 5 b s ) and Rthe are updated to new values. f c is a scale 
factor which reflects the influence of the common error on the fitting. R^QQjys is derived from 



pQCD [2T]. 

The free parameters in the fit are M r , I^ e , 5 r in Eq. ([2]), r r in Eq.©, and Cq, C\, C3 in 
Eq. ([9]). Only relative values of the phases can be extracted, so for simplicity, the phase of 
-0(3770) is set to zero. 



3 Results and discussion 

The resonant parameters of the high mass charmonia determined in this work, together with 
those in PDG2004, PDG2006 and the results given in Ref. [13] are listed in Table [Q The 
updated R values between 3.7 and 5.0 GeV and the fit curves are illustrated in Fig. [TJ The fit 
yields x 2 /d-o.f = 1.06, which indicates a reasonable fit. 

In order to understand the uncertainty of the model-dependence, alternative choices and 
combinations of Breit-Wigner forms, energy dependence of the full width predicted by quantum 
mechanics model |18] or the effective interaction theory |22] , and continuum charm background 
described by a second order polynomial or the phenomenological form used by DASP [6j are 
used. We find that the results are sensitive to the form of the energy dependent total width, but 
not sensitive to the form of background. The DASP background function has six continuum 
production channels, while the effective interaction theory predicts different energy-dependence 
of the hadronic width for different decay channels. However in both cases the best fits give 
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unreasonable results. This may be understood as the inclusive data can not supply enough 
information to determine the relative strength of different decay channels. To understand 
the detailed structures and components of the high mass charmonia states, it is necessary to 
collect data at each peak with sufficiently high statistics, and to develop more reliable physical 
models. This is one of the physics tasks of a tau charm factory, and might be further studied 
with BESIII that is under construction. 

It is worth noting that the change of the resonance parameters affects the effective initial 
state radiative correction factors, and thus affects the R values. Fig. shows comparison 
between the R values published in [14J and the updated values in this work; the difference 
varies with the resonant structure. In general the relative difference is within 3%, and for a 
few energy points the maximum difference is about 6%. Our results are in agreement with 
the previous experiments in most cases, but large differences are observed in some of the 
parameters, such as the mass of the ■0(4160). This is mainly due to the reconsideration of the 
radiative correction factors, and the inclusion of interferences between the resonances. 

Table 1: The resonance parameters of the high mass charmonia in this work together 
with the values in PDG2004 [TT], PDG2006 [12] and K. Seth's evaluations [13] based on 
Crystal Ball and BES data. The total width T tot = T r (M) in Eq.©. 







0(3770) 


0(4040) 


004160) 


0(4415) 


M 

(MeV/c 2 ) 


PDG2004 
PDG2006 
CB (Seth) 
BES (Seth) 
BES (this work) 


3769.9±2.5 
3771.1±2.4 

3771.4±1.8 


4040±10 
4039±1.0 
4037±2 
4040±1 
4038.5±4.6 


4159±20 
4153±3 
4151±4 
4155±5 
4191.6±6.0 


4415±6 
4421±4 
4425±6 
4455±6 
4415.2±7.5 


(MeV) 


PDG2004 
PDG2006 
CB (Seth) 
BES (Seth) 
BES (this work) 


23.6±2.7 
23.0±2.7 

25.4±6.5 


52±10 
80±10 
85±10 
89±6 
81.2±14.4 


78±20 
103±8 
107±10 
107±16 
72.7±15.1 


43±15 
62±20 
119±16 
118±35 
73.3±21.2 


r 

(keV) 


PDG2004 
PDG2006 
CB (Seth) 
BES (Seth) 
BES (this work) 


0.26±0.04 
0.24±0.03 

0.18±0.04 


0.75±0.15 
0.86±0.08 
0.88±0.11 
0.91±0.13 
0.81±0.20 


0.77±0.23 
0.83±0.07 
0.83±0.08 
0.84±0.13 
0.50±0.27 


0.47±0.10 
0.58±0.07 
0.72±0.11 
0.64±0.23 
0.37±0.14 


5 (degree) 


BES (this work) 





133±68 


301±61 


246±86 
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-0.5 - 



-1.0 



_i_ 



X7d.O.f=1.05 




4.7 4.8 4.9 5.0 



3.7 3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 4.6 

Ecm(GeV) 

Figure 1: The fit to the R values for the high mass charmonia structure. The dots 
with error bars are the updated R values. The solid curve shows the best fit, and the 
other curves show the contributions from each resonance Rbwi the interference Ri n t, the 
summation of the four resonances R res = Rbw + Rint, and the continuum background 
R con respectively. 



0.08 






3.7 3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 4.9 5.0 37 3 8 39 40 41 42 4.3 4.4 4.5 4.6 4.7 4.8 4.9 5.0 

Ecm(GeV) Ecm(GeV) 

Figure 2: (I) The comparison of R values between the values published in Ref. [H] 
(triangles: R id) and the updated values in this work (points: R n0 w)- (II) The relative 
differences between the two sets of R values. 
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